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Abstract 

In this paper we classify all 4+1 cosmological models where the spatial 
hypersurfaces are connected and simply connected homogeneous Rieman- 
nian manifolds. These models come in two categories, multiply transitive 
and simply transitive models. There are in all five different multiply tran- 
sitive models which cannot be considered as a special case of a simply 
transitive model. The classification of simply transitive models, relies 
heavily upon the classification of the four dimensional (real) Lie algebras. 
For the orthogonal case, we derive all the equations of motion and give 
some examples of exact solutions. Also the problem of how these mod- 
els can be compactified in context with the Kaluza-Klein mechanism, is 
addressed. 



1 Introduction 

The idea that our world has more dimensions than three is actually older than 
the theory of general relativity itself. Already in 1914, G. Nordstrom^] put 
forward a five-dimensional scalar-tensor theory in an effort to unify gravity and 
electromagnetism. Since it was based upon his own theory of gravitation which 
was soon superseded by Einstein's theory, this work was neglected for several 
decades. 

However, in 1919, T. Kaluza Q constructed a similar unified theory of grav- 
ity and electromagnetism based on the linearized version of the general theory of 
relativity. Kaluza's work, which was published in 1921 and was followed by two 
important papers by Klein |}| Qj, had a very interesting result: five-dimensional 

* S.Hervik@damtp.cam.ac.uk 

1 English translations of this and other related works are found in the book |]J . 
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Einstein gravity could be seen upon as Einstein gravity in four dimensions plus 
electromagnetism. 

The original idea had only one extra dimension, but string theorists today 
believe that we can have up to seven extra dimensions^). Supposedly, six of 
these have to be curled up in a small Calabi-Yau manifold Hence, among 
theorists the question is not if we have extra dimensions, but rather 11 How many 
extra dimensions are there?" and "What is their nature?" 

In this paper we will consider the simplest models, models with one ex- 
tra dimension. We will investigate cosmological models of 4+1 dimensions. In 
3+1 dimensions the classification relies on Bianchi's classification of the ho- 
mogeneous three-manifolds @. Bianchi's work was later generalised to the 
four-dimensional manifolds by Fubini In three dimensions the homoge- 
neous manifolds have a special role. According to the Thurston conjecture, the 
homogeneous manifolds of dimension three are intimately related to the clas- 
sification of three-manifolds fi |l(J. There does not exist a similar conjecture 
in four dimensions; on the contrary, topology in four dimensions is completely 
different. Notwithstanding, we will assume that our 4+1 cosmological model is 
spatially homogeneous. This assumption heavily restricts the number of possi- 
ble four-manifolds to only a finite and manageable number. We know that our 
universe is homogeneous on scales larger than a billion light-years, and thus the 
assumption of homogeneity is by no means a radical and unrealistic one^|. 

The purpose of our study is to investigate the impact of extra dimensions 
on the cosmic evolution of our universe. This may give us some understanding 
of how extra dimensions generically influence the observed 3+1 dimensional 
universe. 

The simplest extension we can think of, going from three to four spatial 
dimensions, is just to assume that we have have a product space: 

£ = M x S 1 (1) 

where £ is a four dimensional spatial manifold and M a three dimensional 
manifold. However, this naive assumption does not necessarily need to be true. 
Experimental data from particle accelerators indicate that the size of this small 
extra dimension must be less than 10~ 18 m. Hence, since the size of the base- 
space M is greater than billions of light-years, the extra dimension may also be 
"twisted" . This twisting is at a global scale, and will therefore be unmeasurable 
in particle accelerators. However, we know that once in the past, the size of the 
universe might have been comparable with the small extra dimensions. Actually, 
it is one of our principal aims to try to explain why one (or more) dimension (s) 
is so incredibly much smaller than the three large ones we see today. 

Some work along similar lines already exist in the literature. Forgacs and 
Horvath jll], [l2| investigated already in 1979 how higher dimensional mod- 
els could influence the cosmology of isotropic FRW universes. This work was 
followed by Chodos and Detweilerjl3| who studied the translational invari- 
ant higher-dimensional cosmological models - also over 20 years ago. Later, 
other people have investigated other models @, [0| ^ 0, Q ^ ® ^ 



22l E3|, B4l E5L E6l E7l |28| . Also, Lorenz-Petzold produced a string of papers 



2 From CMB measurements we know that the fluctuations away from homogeneity are no 
larger than 10~ 5 at the last scattering surface. 
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[§9|, |o[ [HJ |32[ |33j, ||, ||, Hi by lifting the homogeneous models in 3+1 di- 
mensions, one by one, to higher dimensions. Maybe the closest work related 
to this (at least to the authors knowledge) are two articles in the mid-eighties 
p7[ |3a ]. However, there does not seem to be any work which tries to classify all 
the spatially homogeneous spacetimes in 4+1 dimensions and gives a systematic 
approach to the equations of motion. The aim of this paper is exactly to do 
this; give a classification and with the aid of the orthonormal frame formalism, 
we will derive all the equations of motion. 

The paper is organised as follows. First we classify all the multiply transitive 
spaces. We will throughout our paper assume, unless stated otherwise, that our 
space is connected and simply connected. It is of special interest to find all the 
homogeneous spaces which cannot be considered as simply transitive spaces. 
The simply transitive spaces are classified in section ||. This classification heavily 
relies upon the classification of the four-dimensional Lie algebras. We write 
down all the equations of motion and go on and provide with some examples 
and give some exact solutions. Lastly, we go to the question of compactification 
which is one of the key ingredients of the Kaluza-Klein mechanism. 



2 Multiply transitive models 

Let us first classify all the multiply transitive models which are connected and 
simply connected. In dimension three there is only one multiply transitive 
model which cannot be considered as a special case of a simply transitive one. 
This is the well-known Kantowski-Sachs model. The symmetry group in this 
case is four-dimensional and it has three-dimensional subgroup which acts on a 
two- sphere S 2 . 

Our analysis is based on the classification of the homogeneous Rlemannian 
spaces of dimension four due to Ishihara Large parts of our results in 

this section can be extracted directly from this paper. In four dimensions there 
are in all 5 different models which cannot be seen as special cases of a simply 
transitive space. In the following, these will be emphasized.^ 



2.1 Maximally symmetric: dimlsom(S) = 10 

There are three (orientable) maximally symmetric spaces of dimension 4. They 
are the well-known cases: Sr., E 4 and H 4 . These have an isotropy group iso- 
morphic to 5 l O(4). They correspond to the three different FRW cosmological 
models in 5D. Both the Euclidean space and the hyperbolic space have sub- 
groups of the isometry group which acts simply transitive on their respective 
spaces. on the other hand, has no proper subgroup that acts transitively on 
the space at all. Hence, cannot be considered as a special case of a simply 
transitive space. 

3 The homogeneous Riemannian spaces are treated in a more general way in for example 

|40| . However, dimension four seems to be a special case. This is mainly due to the non- 
simplcness of 50(4). 

4 The underlined spaces will correspond to these multiply transitive models. However, they 
will only be underlined in the section were they have the smallest transitive symmetry group. 
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2.2 Kahler manifolds: dimlsom(S) = 8 



Interestingly, there also exist four-dimensional Riemannian spaces which have 
an 8-dimensional isometry groupn. These are the three Kahler manifolds of 
constant curvature: CP 2 , C 2 andfjH 2 . These spaces have an isotropy subgroup 
isomorphic to U(2). Both C 2 and H 2 , have a simply transitive subgroup. 

The Riemann curvature tensor for these spaces can be written in complex 
coordinates as 

R a ps*j = K \9otp95*) + 9aj9sp) ( 2 ) 

where K > 0, K = and K < for CP 2 , C 2 and respectively. Thus these 
spaces have constant holomorphic sectional curvature. The isotropy group is 
the group of motions that leaves the Kahlerian structure invariant. 

2.3 dimlsom(E) = 7 

Spaces that possess a 7-dimensional isometry group come in three classes. The 
first class consists of maximally symmetric three-manifolds times a line: E 1 x 
S 3 , E 4 and E 1 x H 3 . These have an isotropy group isomorphic to SO(3). The 
second category consists of H with a certain 7-dimensional symmetry group. 
The third is the Kahler manifold C 2 with isotropy group SU(2). All of these 
spaces can be considered as a special case of a simply transitive space. 

2.4 dim Isom(S) = 6 

All of these spaces are products of two two-dimensional maximally symmetric 
spaces: S 2 x S 2 , S 2 x E 2 , S 2 x H 2 , H 2 x E 2 , H 2 x H 2 and E 4 . These have the 
isotropy group SO (2) x SO(2). Only S 2 x S 2 does not have a lower dimen- 
sional subgroup of the isometry group which acts transitively on the spatial 
hypersurface. For this model, the metric can be written as 

ds 2 = -dt 2 + a{t) 2 (dd 2 + sin 2 6d4> 2 ) + b{t) 2 (dijj 2 + sin 2 ?/< 2 ) ■ (3) 

2.5 dimlsom(S) = 5 

There are several spaces possessing a 5-dimensional isometry group. Only two of 
them cannot be considered as a special case of a simply transitive space: S 2 x E 2 
and S 2 x H 2 . The isotropy group for both these spaces is 5*0(2). Examples of 
such models for the case S 2 x E 2 case are 

ds 2 = -dt 2 + a(t) 2 (d9 2 + sin 2 Odcf) + b{tf (e~ 2 ^dx 2 + e 2 ^dy 2 ^ , (4) 

and in the case S 2 x H 2 

ds 2 = -dt 2 + a(tf (de 2 + sin 2 9d(f> 2 ) + b(t) 2 (e^'W + e 2 ^+ 2Kx dy 2 ^j . (5) 
The functions a(t), b(t) and j3(t) will be determined by the field equations. 

5 In general for a homogeneous Riemannian space M of dimension n ^ 4 there does not 
exist a closed subgroup of Isom(M ) of dimension r such that ra(n + l)/2 > r > 1 + n(n — l)/2. 
These are the counterexamples for n = 4. 

6 For those who are unfamiliar with the complex hyperbolic spaces, consult for example 
Goldman's book mU. 
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3 Simply transitive models 



In the simply transitive case, we can systematically construct the homogeneous 
spaces using the classification of the 4-dimensional Lie algebras. For a given 
simply transitive space, the Killing vectors obey a certain commutator relation 

= (6) 

where the structure constants are functions of t only. In the following we will 
follow Ellis and MacCallum |Q closely. However we will use a slightly different 
notation: Greek indices (//, v, ...) have range 0-4 over the full five dimensional 
space-time; lower case Latin indices (a,b,i,j, ...) have range 1-4 over the four- 
dimensional spatial hypersurfaces; upper case Latin indices (A, B, ...) have range 
1-3 over three spatially directions. We will assume that the model is orthogonal, 
i.e. we assume that our spacetime is of the form 

S t x M (7) 

where the fluid four- velocity u is orthogonal to St. 

We can now introduce a left-invariant spatial frame e a , which together with 
the fluid four- velocity u = e t forms a orthonormal frame denoted by e^. These 
commute, by definition, with the Killing vectors: 

&,e„]=0. (8) 

The fiinfbein will now obey the commutation relations 

[e^e v ] = 7 P ^ e P- ( 9 ) 

We can relate these commutator function to the connection coefficients for the 
particular orthonormal basis. For an orthonormal basis we define the rotation 
forms by 

de M = e„ ® n" M . (10) 

Denoting the dual frame of by w M , the rotation forms can be related to the 
connection coefficients via 

n% = F>;\ (11) 

Also, they possess the antisymmetry Q^y = — Vt ull . Hence, the connection 
coefficients for the orthonormal frame possess the antisymmetry 



r pfif F^p^ (1^) 



and can be written as 



7 p = _ (yp - r p ) 

I \XV \ pV Vp) 

^ ppv = — — {"fppu + Ipvp ~ lupp) ■ (13) 

The orthogonality condition requires that 

7°0a = 7° a 6 = 0. (W) 
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We can split the remaining part of the structure constants into 

7%, = -0 a b + V a b 

fjk = C%. (15) 

Here are Q ab the angular velocity in the a6-plane of a Fermi-propagated axis 
with respect to the triad e a , and 9 a b the volume expansion tensor. If u M is the 
time-like vector-field orthogonal to the hypersurfaces , then 8 a b is defined by 
d^u — u KV . One can readily see that 9 pu is symmetric and 9 pv u^ = 0. 

We can further split the volume expansion tensor into a trace and trace-free 
part 

0pu = 2 h »v8 + a »v ( 16 ) 

The tensor a pv is the shear tensor and h pv = g pv — u^Uy is the projection tensor 
onto the hypersurfaces S t . 

The structure constants Cy defines the specific Lie algebra under consid- 
eration. These can be separated into a trace (or vector) part and a trace-free 
part 

C % = D % + 5 \ a i - (17) 

where 

ai = ~TV(C*,) ee 

D % = C %-\ Cl W\t% ( 18 ) 

The Lie algebra given by the structure constants is isomorphic to the Lie 
algebra C\j defined in eq. 
The Jacobi identity, 

[e M , [e„, e p \\ + [e v , [e p , e M ]] + [e^, [e M , e„]] = (19) 

reduces now to the following constraint equations 

a k D% = (20) 

D d ab D e cd + D d ca D\ d + D d bc D\ d + D\ b a c + D% a a b + D\ c a a = (21) 
and evolution equations 

d b + X -9a h - (n d b - a d b ) a d = (22) 

D a bc + \eD\ c + 2(n d [b -a d [b )D a c]d + {n a d -a a d )D d bc = 0. (23) 

Equation ( |2l| ) is just the Jacobi identity for C\j in disguise, while eq. (|2fj| ) is 
the trace of eq. (|2l| ) . 
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3.1 The Field Equations 

Let us assume that the energy-momentum tensor is of the form 

Tf,, v = {p + pju^u,, + pg^ + ir^ (24) 

where 7r M „ is the anisotropic stress tensor. This tensor is symmetric and has the 
properties 

7r"„ = = 0. (25) 

The identity = leads to the energy conservation equation 

p+{p + p)9 + n^a^ = 0. (26) 

Together with an equation of state for the fluid, these equations govern the 
evolution of the fluid in our model. 

The Riemann curvature tensor is given by 

R%s - e 7 (r%) - e s (r^ 7 ) + r<\ 7 r% 4 - r« A ,r^ 7 - r^V (27) 

From this expression we can readily calculate the field equations. The 5D Ray- 
chaudhuri's equation is 

()+\6 2 + a ab a ab + ^(p + 2p- A) =0 (28) 

while the (0, a) equations are 

4<r ac a c - a\D\ a = 0. (29) 

The trace of the field equations yields the 5D Friedmann equation (or the 
constraint equation): 

le 2 = \a ab cr ab + Qa 2 + \&+p + k (30) 

3 2 = D\ c D a bc + 2D\ c D b a c . 

The shear equations can be derived from the trace-free part of the field equa- 
tions: 

3 1 
Oab + 0a ab + - (D abd + D bad ) a d + 2a c {a fl b)c + b ab - -h ab b c c = n ab (31) 

where we have set 

b ab = \ (2D c ad D cb d D acd D b cd + 2D\ d D d bc ) . (32) 

Note that 4b a a = — D 2 . The tensor b ab can be interpreted as part of the spatial 
Ricci tensor: 

{i) R ab = b ab + | (D abd + D bad ) a d - 3h ab a 2 . (33) 
Hence, the curvature of the spatial four-surfaces is 

^R = b a a -12a 2 = -^D 2 -12a 2 . (34) 
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3.2 The four-dimensional Lie algebras 

In a series of papers Mubarakzyanov ji^, Q [l5| gave a classification of real Lie 
algebras up to dimension 5.[| A list of all the algebras is given in (47 , 48 which 
is more accessible for non- Russian readers (see also 49 for the four-dimensional 
case). Here we will be interested in the four-dimensional classification and we 
will be using mostly the notation in j|7). A useful reference is also a report by 
MacCallum |§. 

We will investigate Lie algebras from two different point of views; one from 
a geometric point of view, and the other from an algebraic point of view. The 
two different ways of investigating Lie algebras both have their strengths and 
weaknesses and are useful for different purposes. We will start out from the 
geometrical point of view; dividing them into decomposable and indecomposable 
ones. The geometrical picture is completely different in these two different 
classes. The decomposable ones naturally give rise to product spaces^, while 
the indecomposable ones cannot be written as a topological product. 



3.2.1 Decomposable Lie algebras 

Since the three-dimensional Lie algebras are used quite frequently in the litera- 
ture we will use a notation similar to the Bianchi types when the algebras are 
composed of these. The decomposable four-dimensional Lie algebras are 

A 3 © E, and A 2> i © A 2 ,i (35) 

where A3 is one of the Bianchi type Lie algebras, labelled I-IX, and ^2,1 is the 
only non-trivial two-dimensional Lie algebra. In this notation, the extension of 
the type IX Bianchi type will be denoted IXffiR. ^2,1 can be represented by the 
single non-trivial commutator 

[ei,e 2 ] = ei. (36) 

One can show that this Lie algebra acts simply transitively on H 2 . It has been 
known for a while that the commutators for the Bianchi Lie algebras can be 
written 

C A bc = tBCDn DA + S A B a c - 6 A c a B (37) 

where n AD is a symmetric matrix. This is called the Behr decomposition ]5l| . 
The Jacobi identity reduces to the single relation 

n AD a D = 0. (38) 

The Bianchi types can now be classified in terms of the eigenvalues of the matrix 
n AD @. 



3.2.2 Indecomposable Lie algebras 

The classification of the indecomposable four-dimensional Lie algebras are listed 
111 table 0. 

7 The classification of the complex Liealgebras up to dimension 6 was actually done by 
Sophus Lie himself already in the 1890's [46|. 

8 As Lie groups. 



Name 


Commutator relations 


04 


At,i 


[e 2 ,e4] = ei, [e 3 ,e 4 ] = e 2 





-^4,2 


\ei , 64.1 = vb~\ , [e 2 , G4I = e 2 , 
[e 3 , e 4 l = e 2 + e 3 (p ^ 0) 




-44,3 


[ei,e 4 ] = ei, [e 3 ,e 4 ] = e 2 


a i 


v4 4: 4 


[ei , e 4 l = ei , feo, e 4 l = 61 + e 2 , 
fe 3 , e 4 l = e 2 + e 3 




Apq 

^4,5 


[ei, e 4 ] = ei , [e 2 , e 4 ] = pe 2 , [e 3 , e 4 ] = qe 3 
(pq^O, -1 < g < p < 1) 




^4,6 


[ei,e 4 ] =iJei, [e 2 , e 4 ] = ge 2 - e 3 , 
[e 3 , e 4 ] = e 2 + ge 3 (p^O, g > 0) 


s(p + 2g)^ 


-44,7 


[e 2 ,e 3 ]=ei, [ei,e 4 ] = 2ei, [e 2 ,e4]=e 2 , 
[e 3 , ba] = e 2 + e 3 


3 U « 


^4,8 


[e 2 ,e 3 ] = ei, [e 2 ,e 4 ] = e 2 . [e 3 ,e 4 ] = -e 3 





^4,9 


[e 2 ,e 3 ]=ei, [ei, e 4 ] = (1 + q)ei, [e 2 ,e 4 ]=e 2 , 
[e 3 ,e 4 ]=ge 3 (-K q < 1) 


1(1 + 9)^ 


-44,10 


[e 2 ,e 3 ] = e 1; [e 2 ,e 4 ] = -e 3 , [e 3 ,e 4 ] = e 2 





/I 9 

^4,11 


[e 2 ,e 3 ]=ei, [ei, e 4 ] = 2gei, [e 2 , e 4 ] = ge 2 - e 3 , 
[e 3 ,e 4 ] = e 2 + ge 3 , (q > 0) 




^4,12 


[ei,e 4 ]=ei, [e 2 ,e4]=e 2 , 
[ei,e 3 ] = -e 2 , [e 2 ,e 3 ] = -ei 


in 



Table 1: The indecomposable four-dimensional algebras 



Note that some of the parametric limits yields various other Lie algebras. 
We have 

limp_o^S,2 = IV ®R, lim 9 ^o-4^ 5 = VI h ®R, 
lim^o^e = VII h ®R, lim 9 ^o-4^ 9 = A 4)8 , (39) 
lim q ^ A q 411 = A 4 ,i . 

Unfortunately there does not exist a simple expression, as in the three- 
dimensional case, for the structure constants. This is perhaps the greatest 
obstacle to our analysis. However, in several cases, the general form of the 
structure constants can be worked out. 

Let C\j be a special representation of a Lie algebra A. We define the space 
W(A) as 

W(A) = {C%\C% = (A- 1 )" CLA-'.A™, A e GL(4,M)} . (40) 

Two elements in W(„4) correspond to two isomorphic Lie algebras. It would 
be convenient to have a specific parametrisation of W(A) for all of the four- 
dimensional Lie algebras. Unfortunately, such a parametrisation has not been 
found in the four-dimensional case. The union of W(*4) over all possible Lie 
algebras in n-dimensions is called the variety of n- dimensional Lie algebras^ 
In dimensions two and three, the structure of this variety is known, but in four 
dimensions and higher it is not known in detail. 

9 For the case of complex Lie algebras of low dimensions, see [B2|. 
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However, for specific cases, W(^4) can be found explicitly. For instance, for 
the pure vector type algebra, the structure constants are (see Table |l|) 

C% = 8\5% - S k 3 S\. (41) 

Hence, W^A^'g) is the space of all non-zero vectors in R 4 : 

W(A\ : l) - R 4 \ {0}. (42) 

3.2.3 Non-unimodular Lie algebras 

Let us now investigate the algebraic way of looking at these Lie algebras, follow- 
ing MacCallum [ p0[ . We will first investigate the ones that are non-unimodular. 
These are exactly those who have non-zero trace: ^ 0. Let hi be the unit 
vector parallel to ai . For the sake of simplicity, let us assume that our frame is 
orientated so that = aS\. Then we can decompose the structure constants 
into 

c A = e A 



Bi ^ B 

n ab = \c%^ bk a k (43) 
where n ab is a symmetric matrix. The Jacobi identity implies C 4 ab = and thus 

n a4 _ n 4a = q ^ 

Note that the matrix n ab only contributes to the trace-free part of the structure 
constants. Using a "sloppy" notation, but practical for our purposes, we can 
write the trace-free part as 

U BC — e BCDkn a 

D A Bi = Q A B -a6 A B . (45) 

We have still an 5*0(3) orientation we can use to diagonalise n AB . Thus we can 
assume that n AB = diag(n!, 72.2, n.3) by choosing a suitable frame. The Jacobi 
identity reduces now to 

m (2D\ 4 - a) = n 2 (2D 2 24 - a) = n 3 (2D 3 34 - a) = 
n 2 6 3 2 +«3© 2 3 = n 3 Q\ + n 1 9 3 1 = n 1 Q 2 1 +n 2 6 1 2 = 0. (46) 

Let us first assume that the rank of n ab is 3. This leads to D A Ai — 3a/2 which 
is a contradiction, since D a bc is trace-free. Hence, we can assume that TI3 = 0. 

The classification now reduces to finding the eigenvalues n\ and n 2 and the 
matrix 8^ . This is done in MacCallum ; we will take the simplest example 
in this case (even though the dynamical behaviour of these types may be highly 
complex). Assume that the rank of n ab is two. The further analysis splits in 
two cases, determined by the sign of ri\ni. Assume that n\,n2 > 0. This 
leads to D\ A — a/2 and D 2 24 = a/2 (and hence D 3 34 — —a). Further we get 
3 2 = 3 i = plus the constraint n\®\ + U2& 1 2 = 0. This algebra is type 
A 4A2 , which MacCallum calls N22. 

The other Lie algebra with rank(n ab ) = 2 is found when n\ < < n 2 - This 
corresponds to the decomposable algebra A 2 ,i © j4 2 ,i. MacCallum calls this type 
N20. 
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Hence, both of these algebras, At,i2 and ^2,1 © ^2,1, have (after choosing a 
orientation of frame) the following parameters: 

±g\ ey 

I e 2 3 . (47) 

-a 

Doing this analysis for all of the remaining cases (for rank 1 and 0) we can 
find the parameter space for these models. In the appendix, the general form 
of the matrix Q A B is listed for all of the non-unimodular Lie algebras (using 
a particular choice of gauge). For example, in the case where ri\ = n-i = 0, 
cq. ( f46| ) vanishes identically. Thus it remains to classify the different invariant 
properties of the matrix Q A B - These invariant properties determine the Lie 
algebra type. 

3.2.4 Unimodular algebras 

The unimodular algebras are defined by vanishing trace: dj = 0. For these 
algebras the following theorem holds: 

Theorem (Farnsworth and Kerr) For a four-dimensional unimodular Lie 
algebra there will cither exist a p a such that 

C a bc = O a [bPch 6> a = 0, (48) 

or there exists no such p a , and there exists a non-zero l c such that 

C\ c l c = 0. (49) 

A simple and geometric proof of this is given in MacCallum ]50[ . MacCallum 
calls the algebras obeying ( f48| ) and (^jj), Ul and U3 respectively. 

In the class Ul, we can choose an orientation such that p a — p5 A a . Now the 
remaining part to be classified is the matrix 6^. Since the structure constants 
are trace-free, it follows that Q A A = 0. 

The class U3 splits in two; U3I and U3S. Each of these subclasses has two 
members each. The members of the class U3I, and can both be seen 

as parametric limits of a non-unimodular model. The members of class U3S are 
semi-simple and are the decomposable algebras VIII©R and IXffiR. 



AB 



diag(ni,n 2 ,0), D 



B4 



2 





3.3 Examples and solutions 
3.3.1 The types A 3 © E 

It is interesting to investigate the trivial expansion of the Bianchi types explic- 
itly. These models correspond to the simplest Kaluza-Klcin models. They have 
a single extra dimension, and since this dimension is homogeneous, it can be 
compactified into a circle. 

We need to find all the structure constants possible for these models. There 
are two equivalent ways of doing this. Either we can do it purely algebraically, 
or we can do it from a geometrical point of view. We will choose the latter. 

In principle, the extra dimension can be tilted. Let us choose the spatial 
vierbein to have three vectors spanning the vectorspace A3. Hence, we let three 
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vectors have a Bianchi-type algebra, while the fourth vector is orthogonal to 
A%. Let these three vectors which span A3 be denoted e^, B = 1, 2, 3. We can 
now write 

C D AB = e ABC n CD + S D A a B - 5 D B a A (50) 
with the requirement 

a D n CD = 0. (51) 

We know there exists a vector u which commutes with e B . This vector is linearly 
independent of eg so we can write e4 = \ A e A + A 4 u where A 4 7^ 0. Hence we 
get 

[e A , e 4 ] = C D Ai e D = X B C D AB e D (52) 
Thus the structure constants can be written 

C D A 4 = ^ B C D AB , C A ab = 0. (53) 
The trace of the structure constants can now be calculated 
a B — -C lB - -a B 

a 4 - \C\± = ^X D a D = \ D a D (54) 
and the trace-free part is 

t~\A DA { ?A \ 

D bc = e BCDn +-(() B a c - c a B ) 

D A B4 = \ C e BCD n DA -^X A a B + ^S A B a 4 

D\ 4 = a A . (55) 

We still have an 5*0(3) choice of gauge to fix the orientation of the vectors e^, 
hence, we can choose a orientation where the matrix n CD is diagonal. 



Class A models: Let us turn our attention to the Class A models. The class 
A models are characterised by 

a B = 0. 

The Jacobi equations and the (0, i) field equations lead to a set of constraints 
which must be satisfied. Let us choose a frame where n AB is diagonal. The 
constraint equations can now be satisfied with (or a frame can be chosen) 

o A B= = Qab {A ^ B) 
Aa= 

(T4A = QiA- (56) 

Note that we have not assumed a4A = 0, and in general we have not enough 
gauge freedom to put the shear tensor into diagonal form. The case I©R is 
the only case where we can diagonalise the shear vector completely. In the 
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other models there will still remain some off-diagonal shear components. In the 
absence of anisotropic stress, we will get further constraints from the off-diagonal 
shear equations. These can be satisfied with 

= = (74A- (57) 

Thus, in this case the shear must be diagonal. Henceforth we will assume that 
this is the case. 

The tensor b a b reduces to the form 



1 



b AB = 2n Ac n c B - n c c n AB + ^Sab (n c c ) 2 -2n AC 
b ia = (58) 
The independent shear equations are now (ct 4 4 = — a A A ) 

oab + 0a AB + b AB - ^5 AB b c c = (59) 

while the remaining Jacobi equations are 

n AB + Un AB + a c c n AB - 2a {A c n B ^ c = 0. (60) 

Note that these are only valid in the choice of gauge mentioned above. 

Example: A II©M perfect fluid solution. Let us consider a specific case. 
Assume that the matrix n AB has only one non-zero eigenvector. We assume 
that 

n AB = diag(n, 0, 0) (61) 
which corresponds to the II®M Lie algebra. We introduce new variables o\ 2 3 

by 

an — 3<7i 

C22 = — (Tl — + V6CT3 

033 = -ai- V2a 2 - V60-3 

(T44 = — (Tl + 2V2ct 2 . (62) 
The equations of motion now reduce to 

5 

a 1 + Qoi + —n 2 = 
24 

& 2 + 6<t 2 + -^-j=n 2 = 
6v6 

0"3 + 0(73 = 

n + -On - 5am - 2\/2cr 2 n = 0. (63) 
4 

In addition to these, Raychaudhuri's equation and the energy conservation equa- 
tion must be fulfilled. We will assume in this example that the perfect fluid obeys 
a 7-law equation of state: p = (7 — l)p. 



,C 



n 



B 
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There is a specific case where we can solve the equations of motion exactly. 
We can search for a self-similar solution with the properties that ci ; 2,3 oc £ -1 , 
p oc t~ 2 and 9 oc By doing this, we can find a solution which is the II©K 
version of the Collins-Stewart type II perfect fluid solution. This solution is 
given by 

<ri = --|-(2 7 -l)i- 1 
067 

t 3 = 

e = -r 1 (64) 

7 

where l/2<7<2.0 The metric can written as 

ds 2 = -dt 2 +t^ 8 - 5 ^ (dx+-ydz 



7 

+t T^(3T+4) ( dy 2 + d ^ + t ^( 6 - 7)dw 2 ( 65 ) 

where c 2 = (2 - 7) (27 - 1)/H. Note that if 7 > 8/5 - if the matter is stiff 
enough - one of the dimensions will contract to an arbitrary small size. 

3.3.2 Solutions to the model 

This model is of pure vector type and is particularly easy. We have already 
worked out the space of structure constants for this model. The structure 
constants can be characterised by a non-zero vector ai in R 4 . We choose an 
orientation of the frame so that ei is aligned with a,. Hence, ai = 5\a for 
a ^ 0. 

The Jacobi equations ( f2p|) and (pi] ) are satisfied by construction, and eq. 
23]) is trivially satisfied. Equation (g2^ leads to 

a H — 6a + ana = (66) 
4 

and 

^I2=ci2, 13 =(Ti3 and Q 14 = a u (67) 

Eq. (^9|) equation tells us that 

a al a = 0^a al =0. (68) 

Hence, the shear has only non-zero components in the 2, 3 and 4 directions. We 
still have an unused freedom of choosing the orientation of the vectors for 
i = 2, 3 and 4. Hence, we can choose a frame where the shear a a }, is diagonal. 

°It should be noted that in 4+1 dimensions, radiation has 7 = 5/4, and the strong energy 
condition requires 7 > 1/2. Hence, inflation can occur for 7 < 1/2. 
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Henceforth, we will also assume A = = 7r a {, and that the fluid obeys the 
barotropic equation of state p = (7 — l)p. 

The dynamical systems approach has proved to be a powerful tool in cos- 
mology (see for example |53|). Let us use this method to solve the complete 
system of equations in the model. We parametrise the shear with 

a ah = diag(0, —2a + , a + + V3a_, a + - V3(T_) (69) 

and introduce a new time coordinate with 

dt _ 4 
H 1 ~ 9' 

Introducing the expansion normalised variables 



(70) 



_2V2a± 4a 8p 

S±-— ^— , A- T , il-— (71) 



the equations of motion can be written as 



54 = (?-3)S± 
A' = qA 

n' = 2[g-(27-l)]n (72) 



where 



q = 3S + (2 7 -l)fi 

S = + Y?_ 

1 = E + ^ 2 + r2. (73) 

The last of these is the constraint equation, and q is the deceleration parameter 
defined by & = -(l + q)9. 

The system of equations turns out to be quite simple, which makes it possible 
to solve these equations exactly. The solutions are 



£+ = 



P±e~ 3T 



1/2 



n = 



(42 e -6r + e -2(27-l)T + #2) 
e -2(2 7 -l)r 

A 2 e- 6T + e -2(2 7 -i)r + x 2 

A = (74) 

(A 2 e~ 6T + e -2(27-i)r + K 2 ) 1/2 

where p\ + p 2 _ = A 2 . Note that the solutions are Kasner-like near the initial 
singularity. At late times, the universe approaches the 5D Milne universe if 
1/2 < 7 < 2 and the flat 5D FRW if < 7 < 1/2. 
We introduce the functions 



W(i) s (A'+i'-l + KVY 
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The metric for the solutions can now be written 

-x/2co S (0)/3(t) j 2 



e 



+e V2cos(0+ 7r /3)/3(t) dz 2 + e v / 2cos(0- 7r /3)/3(t) rfM; 2 



(76) 



3.3.3 The four-dimensional Nil 4 case 



The only indecomposable Lie algebra in four dimensions that is nilpotent, is the 
type -44,1. This algebra is the Lie algebra of the four-dimensional nil-geometry, 
denoted by Nil 4 . For the sake of illustration we will determine the possible 
structure constants for this model. 

What we would like to have is a representation of the space W(Ai,i). But 
finding the whole space is by no means necessary. We have a gauge freedom 
which corresponds to an 0(4)-rotation of the spatial frame. By the Gram- 
Schmidt process, we can show that every A € GL(n, R) can be written as 

A = P • R (77) 

where R G 0{n) and P E PT(n) where PT(n) C GL(n,M) is the group of 
all upper-triangular matrices with positive entries along the diagonal. Hence, 
by choosing a suitable frame, we can "gauge away" the rotation matrix R. It 
suffices therefore to look at the space 

V(A) = W(A)/0(A) = {C%\C% = (P" 1 )* CP^P™, P e PT(4)} . (78) 

PT(4) is a Lie group, thus if P E PT(4) then P" 1 E PT(4). From Table we 
have 

C x 24 = l, C 2 34 = l. (79) 

By calculating P(A) we find that there are essentially only three non-zero com- 
mutators. These are 

C\ 4 , C\ 4 and C 2 34 . (80) 

(see also the appendix where the A^^ can be seen as the a — > limit of 4^4.) 
The Lie algebra is trace-free, so C*y = D\j, Note that in the general case, 
the three-dimensional Ricci tensor is not diagonal. Hence, if the shear is to be 
diagonal, then we have to set either D 1 2 4 or D 2 3i to zero if no anisotropic stress 
is present. Let us assume therefore that ir a b = 0, C 34 = 0, fl a b = and a a b 
diagonal. 

We introduce the two curvature variables 

Nt = D\ A , N 2 = D\ A . (81) 

So in this model, there will be 5 variables left to describe the geometry of the 
spatial hypersurfaces. These are 

Ni, N 2 , (71,2,3 (82) 
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where cr 1,2, 3 are defined in cq. (|62j). Together with the matter equations, the 
evolution equation for these variables can now be written down. 

In the case of a 7-law perfect fluid, one can also find a specific self-similar 



solution. The metric for this Nil 4 solution is 



ds 2 = -dt 2 + t^ 4 -^ (dx + ^ydw^j 2 + if ^ (^dy + ^zdw^j 2 

+i^+ 2 ^z 2 +^ +1 W (83) 

where c 2 = (2 - 7) (27 - 1)/10 and 1/2 < 7 < 2. The matter-density of the 
perfect fluid is 



For this solution, if 7 > 4/3 then there will be one contracting direction. How- 
ever, this solution is only a particular solution, corresponding to a fixed point 
in the dynamical system. The more general behaviour and the nature of this 
solution for the l\lil 4 -world will be the subject of a future work. 

3.4 Exceptional models A\ 

For the Bianchi models in 3+1 dimensions, there is one exceptional case VI^ , g 
for which one of the i?oa-equations vanish identically; it can have one extra 
shear degree of freedom. 

In the 4+1 dimensional case, we have many such exceptional cases. By 



inspecting eq. (29), we see that we can have an additional shear degree of 



freedom in the following cases: 



^4,2 : 


P 


API . 
A i,5 ■ 


<1 


API ■ 
^4,6 • 


P 



-1 

1+p 



2 



(or with q,p interchanged) 



We will denote these models which have this additional shear degree of freedom 
with A 



with A\. Note that the case VT_ 1 ^ g © K is obtained in the limit of A p ^ as 



lim Af: = VI* 1/9 8 R. (85) 
The other exceptional models are special for the 4+1 dimensional case and 

-A —A* 

have no 3+1 dimensional analogue. Note also that in the cases A 4 |' 3 and 
A 1 ^ 1 * we can have two additional shear degrees of freedom. However, the total 
number of parameters remains the same because in these cases we have one less 
commutator function (see appendix)^. 

Hence, for a given Lie algebra, these exceptional models have more degrees 
of freedom than any of the other models. 

11 This is because we have an additional symmetry in these cases which makes it possible to 
gauge away one more commutator function. 
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4 Compact ificat ion 



As we now have presented all the homogeneous models, we note that most of 
the models are non-compact. The Kaluza-Klein mechanism needs a small and 
compact dimension to work, and hence, all of the models which do not have 
a compact dimension cannot be a proper model for the Kaluza-Klein mecha- 
nism. However, in many of the models we can construct compact versions of 
the models. In a Kaluza-Klein model one usually have a product space 

M x N (86) 

where M is the 3+1 dimensional base space and N is a compact space. For ex- 
ample, in string theory N is a Calabi-Yau manifold with 3 complex dimensions. 
In our model, which has in total 5 dimensions, ./V must be one-dimensional, and 
thus N = S 1 . 

However, one could instead imagine a more general model where the extra 
dimension is small, compact and perhaps "twisted". Hence, we do not neces- 
sarily need to constrain ourself to a product space, but we can allow for a more 
general space. We can assume that our space is a fiberbundle P with compact 
fibers. For a fiberbundle there will also exist a projection map n : P t-^> M. M is 
called the base space, and in our context this is our four-dimensional spacetime. 
We demand that for a p e P we have 

7T- 1 (p)=S 1 . (87) 

Locally we will not see the difference between these "twisted" spaces and the 
product spaces; for any p <G P there will exist an open neighbourhood U £ M 
such that 

U x S 1 C P. (88) 

The difference will be at a global scale. In a cosmological setting, the space M 
will be of cosmological size, and hence, we have to go to cosmological scales to 
see the difference between the twisted and the un-twisted versions. Notwith- 
standing, in the very early universe the cosmological scales may have been com- 
parable to the small extra dimensions. Thus these twisted spaces may have had 
an important impact on the early evolution of the universe. 

To construct a compact dimension we can proceed as follows. We find a 
discrete subgroup T C Isom(P) which acts freely and properly discontinuously 
on P. We identify now points p and q for which there exist a 7 e T such that 
p = j(q). There are usually many such groups for a homogeneous manifold, but 
not always. Also in many cases, the compactification radius can vary as we move 
along the base space. For the simply transitive models this compactification can 
always be done. These models have four linearly independent Killing vectors, 
£j, which act freely on the spatially homogeneous hypcrsurfaces. Hence, for 
any Killing vector £ which is a linear combination of the ^'s we can - through 
exponentiation <p = cxp £ - find such a freely and properly discontinuous acting 
group T generated by the element <j>. 

In the decomposable models, the compactification depends on the single 
extra dimension, which is trivial. In many of these cases, we can also compactify 
the base space as well, and in three dimensions this has some very interesting 
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consequences [pll |5l |6l |7l 58, 59 . In particular, due to Mostow's rigidity 



theorem, all compact hyperbolic spaces of dimension three or higher, are rigid; 
i.e. they do not allow an anisotropic expansion. For the Bianchi types III and 
VIII (and hence types III©K and VIII+R as well) the effect of compactification 
is quite the opposite; they can have a unbounded number of free parameters. 
The compactification induces so-called moduli parameters which will increase 
the number of free parameters. For example, the vacuum type I©R model has 2 
free parameters, while the compactified version of type IffiK into a four-torus T 4 
has 2+16 free parameters^ . This may have a significant effect on a quantum 
theory of gravity and on quantum cosmological models. 

We would leave the compactification question open, but we will provide with 
an example in which we compactify both the base space and its fibers. 



Compactification of the nilpotent geometries. We start by considering 
a three-dimensional torus T 3 . The torus can be constructed from the usual E 3 
by identification under a discrete group of translations. 

Consider the product between the torus and the finite interval [0, 1] 

T 3 x[0,l]. (89) 

The torus can be viewed upon as the unit cube in E 3 with the usual identification 
of the boundary. We will now identify the two tori on boundary of (^) as 
follows. Take a matrix A € SL(3, Z) (the special linear group with integer 
entries). <SX(3,Z) is the mapping class group of T 3 , hence, the mapping A(T 3 ) 
maps the torus isometrically onto itself. To obtain a compact space we can 
therefore identify the boundary T 3 x {0} with T 3 x {1} under the mapping of 
A. Thus if p, q G T 3 then 

{?} x {0} ~ {<?} x {1}, iff A(q) =p. (90) 

We have now constructed a compact manifold M = T 3 x [0, 1]/ ~. If the matrix 
A is the identity matrix, then M = T 4 and hence, is of type / ®R. The three 
other possibilities are the two nilpotent spaces^ Nil 3 x S 1 and Nil 4 , and the 
solvable case Sol 3 x S 1 . 

We define the characteristic polynomial p&( A) by 

Pa (A) = det (A - Al) . (91) 

The roots of pa(A) determine whether we have a nilpotent group or not. 

We have a nilpotent geometry if all three roots of A are equal to 1. If this 
is not the case, then we have the solvable case Sol 3 x S 1 . All of the nilpo- 
tent geometries in four dimensions can be compactified completely in this way, 
depending on the reducibility of the matrix A. 



5 Conclusion and Outlook 

In this paper we have classified all spatially homogeneous cosmological models 
of dimension 4+1 where the spatial hypersurfaces are connected and simply 

12 For the 3+1 dimensional case this is illustrated in jsc|] - 

13 Nil 3 X S 1 has the type II ffiM while the Sol 3 X S 1 has a type VIo © R Lie algebra. 
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connected. We found five multiply transitive models which cannot be seen upon 
as special cases of the simply transitive models. These are in some sense the 
four-dimensional versions of the Kantowski-Sachs (KS) model. Three of these 
models were directly linked to the KS case since they consisted of products of 
manifolds where at least one of the components was a sphere S 2 . 

For the simply transitive models, the classification of the four dimensional 
Lie algebras provides us with all the possible models. Among these models we 
should expect many interesting phenomena. For example, they may give us 
an understanding of how extra dimensions can affect the evolution of our four 
dimensional universe. The idea of extra dimensions is by no means a new idea, 
but higher-dimensional cosmologies seems to be little understood. 

A special class of solutions to the field equations are particularly interesting. 
The plane-wave solutions are known to be solutions describing gravitational 
waves propagating through spacetime. A feature of these solutions is that they 
possess an extra symmetry in addition to those arising from the requirement of 
spatial homogeneity. These solutions have a null Killing vector. The total sym- 
metry groups are therefore higher than the generic solution of the homogeneous 
field equations. These solutions also have a particular interest in string theory 
because they admit supersymmetry. In 4+1 dimensions there exists a large class 
of plane wave solutions. As a matter of fact, there is a five parameter family of 
vacuum plane- wave solution in 4+1 dimensions. Some of them generalises the 
known plane wave solutions in 3+1 dimensions, others are new and special for 
4+1 dimensions. These solutions will be the subject of a future work. 

A recent work which should be mentioned in relation to this, is a work 
by De Smet (6lJ. In this work, a Petrov classification of algebraically special 
five-dimensional spacetimes is constructed. However, what special role these 
spacetimes may have for multidimensional cosmology is not known. 
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Appendix: The structure constants for the inde- 
composable Lie algebras 

We will here give a list of all the structure constants for the indecomposable 
ones. We have chosen a specific gauge (orientation of frame) to simplify the 
expressions and clarify their structure. We write them in terms of the two 
matrices n AB = diag(ni, 712, 0) and C A Bi = Q A b- Note also that the obtained 
bound for the structure constants may in some cases not hold globally. 

rank(n' 4s ) = 2 : m, ri2 ^ Q 



" 3 

2 




B — 



(92) 
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• ^4,12 : o < m,Jl2 

• A 2 ,i ®A 2>1 : m < < n 2 . 
rank(n AS ) = 1 : m 7^ 0, ?i 2 = 



e- 



| a e 1 , eV 
e 2 2 e 2 3 
e 3 2 e 3 3 



. ^ )9 : e 3 2 = 0, e 2 2 > 0, ge 2 2 = e 3 3 , e 2 2 + e 3 3 = fa 



q = 1 =>- <d A 



0. 



• v4 4 , 7 : e 3 2 = 0, e 2 3 > 0, e 2 2 = e 3 3 > o 

• A\ n : Tr(T) = |a, 4g 2 det(T) = (Tr(T)) 2 (1 + g 2 ) where 



T = 



e 2 2 e 2 3 
e 3 2 e 3 3 



(93) 



(94) 



rank(n AB ) = : ti\ = n 2 = 



e\ e 1 , 
e 2 2 e 2 3 
e 3 o 



& A A = 3a. 



e 1 



(95) 



. Al%-. e\ > o, P e\ = e 2 2 , g e\ = e 3 3 . 
p = 1 =>. q\ = 0, P = 9 =► e 2 3 = 0, P = q = 1 

. ^ )2 : e 2 2 = e 3 3 f o, e 1 ! = P e 2 2 , e 2 3 > o, 

• a 4 , 3 : e 2 2 = e 3 3 = 0, e 2 3 > o 

• A iA : Q\ = 6 2 2 = e 3 3 = a, 9 2 3 > 

• A p 4 q 6 : e\ = 3ap/(p + 2q) £ 0, Tr(T) = 6aq/(p+2q) > 0, 4g 2 det(T) = 
(Tr(T)) 2 (1 + g 2 ) where T is defined in eq. @. 

The rest of the non-unimodular Lie algebras can be obtained by taking various 
parametric limits of the above algebras. 
Letting the trace go to zero, leads to 



A 



4,4 





-44,i 



(96) 



and hence, all the indecomposable Lie algebras and the non-unimodular algebras 
can be extracted from the above. The remaining ones, the decomposable of class 



A, are included in the analysis in section 3.3 
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